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摘要 : 磁 流体 数值 模拟 是 空间 物理 研究 的 重要 手段 。 本 文采 用 具有 TVD (Total 
Variation Diminishing ) 特性 的 Lax-Friderichs 差分 格式 求解 了 GLM-MHD 


(Generalized Lagrange Multiplier- Magnetohydrodynamics) 方程 。 为 降低 格式 的 
数值 耗 散 ， 引 入 耗 散 修 正 系 数 对 算法 的 通 量 计 算 过 程 进行 了 改进 。 对 二 维 Rotor 
算 例 和 磁 云 -电流 片 相互 作用 算 例 的 模拟 结果 表明 , GLM-MHD 方法 可 以 有 效 的 控 
制 磁场 散 度 误差 相对 于 泊 松 校正 法 可 以 节省 一 半 以 上 的 计算 时 间 。 在 不 破坏 格 
式 稳定 性 的 基础 上 ， 耗 散 修正 系数 有 效 降 低 了 算法 的 数值 耗 散 。 
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Improvement and application of Lax-Friderichs scheme in 


MHD numerical simulation 

Abstract: Magnetohydrodynamics (MHD) numerical simulation is an important tool 
for space physics research. In this paper, Lax-Friderchs scheme with TVD property is 
employed to solve GLM-MHD equations. The diffusion turning coefficient is 
introduced for scheme optimization. Simulation result of 2D rotor test and magnetic 
cloud current sheet interaction test demonstrates GLM-MHD method's divergence 
control capability. The simulation consumes less than half of the computational time 
comparing with simulation utilizing Poisson correction method. While numerical 
stability is not damaged, numerical diffusion is reduced by the diffusion tuning 
coefficient. 
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日 地 空间 物理 学 中 的 许多 现象 可 以 用 磁 流 体力 学 (Magnetohydrodynamic， 
简称 MHD) 模型 来 描述 。 由 于 MHD 方程 组 的 高 度 非 线性 ， 很 难 求 出 一 般 问 题 的 
解析 解 。 应 用 MHD 方程 组 进行 理论 研究 时 , 一 般 要 借助 数值 模拟 的 方法 。 同时 ， 
MHD 数值 模拟 可 以 给 出 空间 物理 学 现象 在 所 有 关心 的 空间 位 置 和 时 间 点 上 的 物 
理 参 数 ， 弥 补 了 观测 数据 在 时 间 与 空间 上 的 局 限 。 此 外 , 空间 物理 学 所 研究 的 现 
象 ， 很 难 在 实验 室 中 重复 出 来 。MHD 数值 模拟 为 空间 物理 理论 研究 提供 了 可 控 
的 实验 途径 。 随 着 数值 技术 和 计算 技术 的 不 断 发 展 ，MHD 数值 模拟 在 日 蝎 、 行 
星际 、 磁 层 以 及 电离 层 等 大 尺度 .复杂 物理 过 程 研 究 中 得 到 了 越 来 越 广泛 的 应 用 。 
[1] 

MHD 方程 是 Navier-Staokes 7; f£. Maxwell 方程 以 及 广义 欧姆 定律 的 耦合 ， 
具有 典型 的 非 线性 以 及 非 严格 的 双 曲 性 。MHD 方程 的 数值 模拟 方法 ， 大 都 从 流 
体力 学 (Hydrodynamic, 简 称 HD). 方程 的 数值 模拟 中 借鉴 而 来 。 然 而 ，MHD 方程 
的 特征 系统 比 HD 方程 的 要 复杂 很 多 [2]， 这 就 使 得 数值 求解 MHD 方程 组 比 求解 
HD 方程 组 要 困难 。 早 期 数值 求解 MHD 方程 的 方法 主要 是 以 Lax-Wendroff (LF) 
格式 为 代表 的 有 限 差分 方法 。[3] 虽 然 他 们 的 精度 可 以 达到 三 阶 ， 但 他 们 在 间断 
附近 会 产生 振荡 ,不 能 精确 的 捕捉 到 激 波 、 接 触 间断 等 非 线 性 结构 .1983 年 Harten 
首次 提出 了 数值 格式 的 总 差 不 增 (Total Variation Diminishing， 简 称 TVD) 性 质 
(Harten1983)[4]。 有 共 备 这 种 性 质 的 数值 格式 能 够 保证 数值 解 的 单调 性 ， 有 效 地 抑 
制 间断 附近 的 振荡 ， 一 般 被 称 为 TVD 格式 。 早 期 的 TVD 格式 在 极 值 点 上 只 有 一 
阶 精度 ， 通 过 MUSCL (Monotonic Upstream-Centered Scheme) 重 构 方 法 ， 可 以 
使 格式 的 精度 达到 二 阶 [5]。 通 过 PPM (Piecewise Parabolic Method) 重 构 方法 ， 
可 以 使 格式 的 精度 达到 三 阶 。[6] ENO/WENO Weighted Essentially Non-oscillatory ) 
重 构 方法 [7]， 能 使 格式 的 精度 达到 四 阶 乃 至 更 高 。 然 而 ， 由 于 高 阶 重 构 方法 在 
间断 附近 会 自动 降低 精度 以 抑制 谋 荡 , 对 于 存在 大 量 非 线性 间断 的 问题 , 没有 必 
要 使 用 精度 在 二 阶 以 上 的 重 构 方法 。[8] Tanaka 等 提出 了 求解 MHD 方程 的 一 种 
TVD 格式 ， 并 将 磁场 分 裂 成 本 底 磁 场 场 Bo) 和 变化 磁场 (B1)， 以 增强 格式 在 
地 球 磁 层 模 拟 中 的 性 能 。[9] 汉 学 尚 、 魏 泰 思 以 及 范 全 林 等 采用 修正 的 
TVD-Lax-Friedrich 格式 ， 发 展 了 一 种 快捷 共有 TVD 特性 的 组 合 数值 方法 [10-13]， 


在 三 维 球 坐标 以 及 2.5 维 直 角 坐 标 系 下 对 MHD 进行 模拟 ， 解 释 并 证 明了 大 尺度 


太阳 风 结 构 以 及 磁 云 边界 等 问题 。 


Maxwell 方程 组 要 求 磁场 的 散 度 必须 处 处 为 零 (V .8 = 0)。 但 在 多 维 的 MHD 


数值 模拟 中 ， 这 个 条 件 的 满足 需要 特别 的 处 


理 方法 。 这 是 MHD 数值 模拟 比 HD 


数值 模拟 更 复杂 的 男 一 个 原因 。 在 保持 磁场 散 度 为 0 的 方法 中 , 常用 的 泊 松 校正 
法 [13] 需 要 在 每 一 步 求解 泊 松 方程 ， 数 值 计算 耗 时 比较 严重 。8 波 法 [14,15] 则 会 


破坏 MHD 方程 的 守恒 性 。 CT CConstrain Transport) 一 般 需要 将 交 


普 网 格 来 实现 ， 


步骤 比较 繁琐 。[6] 广 义 拉 格 朗 日 滋 子 (Generalized Lagrangian Multiplier， 简 称 
GLM) 方法 [16] 通过 在 MHD 方程 中 增加 变量 y 的 方式 , 构成 GLM-MHD 方程 组 ， 


在 不 破坏 方程 守恒 性 的 前 提 下 , 使 磁场 散 度 误差 名 


E 输 运 的 同时 被 耗 散 。 标 准 算 例 


的 测试 结果 表明 ,这 种 方法 能 够 有 效 的 控制 磁场 散 度 误差 , 得 到 精确 的 数值 结 


同时 可 能 在 间断 附近 产生 不 正确 的 解 。 


本 文 用 到 修正 的 TVD-Lax-Friedrich 格式 求解 了 GLM-MHD 方程 ， 


引入 了 耗 散 


修正 系数 (diffusion reduction coefficient) 以 降低 格式 的 耗 散 。 通 过 对 二 维 Rotor 


问题 和 行星 际 磁 云 边 界 层 问 题 的 模拟 ， 验 说 


二 、 基 本 方程 和 方法 
1、 控 制 方程 


本 文 求解 的 是 2.5 维 守恒 形式 的 电阻 GLM-MHD Z7 f: 


QU , OF 


E 了 格式 的 性 能 。 
TVD-LF、TVD-TD 格式 求解 GLM-MHD 方程 的 具体 步骤 。 第 三 、 
问题 和 行星 际 磁 云 边界 层 问题 的 模拟 结果 。 第 五 章 得 出 结论 。 


本 文 第 二 章 介绍 了 


四 章 分 析 了 Rotor 
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HELO) = 2 775 为 拉 普 拉 斯 (Laplacian) 算 符 。 方程 中 的 拉 格 朗 日 乘 子 y 是 
一 个 标量 ， 它 是 与 磁场 散 度 误差 相关 的 一 个 量 。cn 用 来 散 度 误差 传播 速度 。 若 用 
Cfa Cfy 表 示 x 和 yy 方向 的 快 磁 声 波 速 度 , 则 可 令 ch = max(]ul + crx, Iv] + cry) 
从 而 在 不 对 计算 时 间 步 长 产生 附加 限制 的 前 提 下 , 以 最 大 的 速度 将 产生 的 磁场 散 
度 误差 输 运 出 去 。 同 时 ， 按 照 文献 [16]， 将 常数 cp 的 取 值 为 cp = V0.18cn。 
2. GLM-TVD-LF 格式 
我 们 将 冯 学 尚 等 [10] 求 解 MHD 方程 的 TVDLF 格式 应 用 到 求解 GLM-MHD 方程 
上 。 这 种 格式 将 一 个 时 间 步 的 求解 过 程 分 成 预 估 步 和 校正 步 两 个 步 又。 求解 的 具 


体 过 程 为 
CG1) 预 估 步 : 通过 预 估 步 ， 得 到 U 在 两 个 时 间 步 之 间 的 中 间 值 
1 
nt+= At 1 1 At 
U d — UL [F (ut; $ -8Up) -F (Uf. - 28u )| E | (Uf. 4 
1 1 At 
z8Up) — G (Ut, —26U7)| + 5s UD (6) 


(6) PF, G, 的 具体 形式 分 别 见 〈3)，(4) XX. 


SU} = minmod(AUP 4, jy, AU ajo jk) AU aje = Uaj Uik (7) 


SUP 的 定义 与 之 类 似 。 这 里 使 用 minmod 限制 器 来 抑制 间断 附近 震荡 的 产生 ; 


minmod(x, y) = sgn(x)Max(0, min[|x|, y sgn(x)]) 
1, x»0 
sgn(x) 2410, x20 (8) 
—1, x «0 
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以 体现 源 项 Sz 的 作用 


(2) 校正 步 
通过 校正 步 ， 得 到 物理 量 在 下 一 个 时 间 步 的 值 U? 宫 
mi+1 rm _ At an+ Oo pnt _ At an+ u 
"Lk = Ug s 2537: un Ay (e, 
o 
以 x 方向 的 通 SPUR 绍 校正 步 中 通 量 的 计算 方法 。 定 义 
L — rr”+1/2 n R mt+1l/2 1d ep 
i. zik = T3 =U; Ui ik LIE 一 70Ui+1 
其 中 ，6UFP 的 值 已 在 预 估 步 中 由 (7) 式 算出 ， 这 里 无 需 再 次 计算 。 


_ Bxiryjk — 
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在 计算 (6) 式 后 ， 参 考 文献 [16] 中 的 方法 对 炒 值 了 


以 Wi; 7 作为 % 在 校正 步 中 使 用 的 值 。 
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在 计算 Pajx 时 ， 参 考 [16] 的 方法 ， 首 先 解 出 GLM 子 系统 的 黎 曼 问题 ， 得 到 
L 
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pa = ip Vuit en (pR L (11) 
i+} j,k 2 2 \ xitLjk xi+Ż jk 
2 2 2 
m^ pL — pR — L — R y 
HS By iik By aii By ii pg? Vitii Wg Wa 随后 ， 用 
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Lax-Friedrich 方法 计算 出 通 
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enira) (ou) 
定义 crx(U) 为 以 UV 中 物理 量 的 值 计 算出 的 


UF, 
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在 计算 完 (9) 式 后 ， 采 用 与 预 估 步 相似 的 方法 ， 更 新 的 值 


c2 
nAt 
*N 二 1 tle 
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p 
3. GLM-TVD-TD (Tunable Dissipation) 修正 
用 LacFriedrich. 方法 计算 数值 通 量 ， 虽 然 稳定 性 较 好 ， 但 其 将 粘性 项 
lAmaxlır ( UP, 一 UL, 的 系数 固定 为: 数值 粘性 较 大 。Diebel、 Yalim[17,18] 
Zj, 


it; „j,k 
提出 引入 一 个 取 值 范围 在 (0,1] 之 间 的 耗 散 修正 系数 =， 用 来 控制 通 量 的 耗 散 的 大 
小 。 于 是 ， 数 值 通 量 (13) 可 以 改写 : 


1 
Ana 1 L R E R L 

F 2 --|F[(U F| U; 一 =|4 U 一 U 
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的 值 取 的 越 小 ， 数 值 通 量 的 耗 散 越 低 。 但 e 取 的 过 小 时 ， 间 断面 附件 会 
生 数 值 振荡 ， 使 得 计算 不 能 稳定 进行 。 目前， 非 线 性 方程 数值 解法 的 稳定 性 还 没 
有 严格 的 判断 准则 ，gs 难 以 通过 理论 推导 给 出 。 我 们 采用 Yalim[18] 的 方法 ， 通 过 
采用 不 同 s 值 进行 多 次 数值 实验 的 方式 ， 选 取 结 果 中 不 存在 振荡 等 非 物 理 现象 的 
最 小 的 e 值 。 这 样 就 可 以 在 追求 精确 度 和 保证 数值 计算 稳定 性 之 间 寻 求 到 一 个 最 
佳 平衡 点 。 
三 、 二 维 Rotor 问题 数值 试验 

为 了 对 GLM-TVD-LF 和 GLM-TVD-TD 方法 进行 验证 ， 本 文选 取 Rotor 问题 进 
行 数值 模拟 。1999 年 Balsara and Spicer[19] 为 了 研究 强 扭转 Alfven 波 结构 和 传播 
提出 了 二 维 Rotor 问题 。 由 于 Rotor 问题 具有 典型 的 磁场 和 流 场 间断 面 结构 ， 使 
其 成 为 了 验证 MHD 模拟 方法 的 重要 算 例 。 

本 文选 取 的 模拟 区 域 为 x,y € [一 0.5, 0.5]， 采 用 400 x 400 的 均匀 网 格 。 算 例 的 初 
始 条 件 [16] 为 : 


(10, —wy, wx) if r € rg 
(P, v vy) = (1 +9f,— cem em) Ifry«Tr «mr (16) 
(1, 0,0) ifr2n 


E Pw = 20, ro = 0.1, r = 0.115, r = 4x? + y?, f-2(n-r)(m-r). BS 
B, = 5.0/V4r， 热 压 p = 1， 绝 热 系 数 y = 1.4. ift = 0.15 时 停止。 


图 1 中 自 上 而 下 分 别 为 使 用 [10] 中 的 TVD-LF 方法 、 本 文 的 GLM-TVD-LF 方法 
和 GLM-TVD-TD 方法 模拟 得 到 的 Rotor 问题 的 密度 p 等 值 线 和 磁场 B? 等 值 线 图 ， 
其 中 耗 散 修正 系数 选取 为 通过 试验 得 到 的 最 优 值 s = 0.5。 通 过 对 比 等 值 线 结果 ， 
我 们 发 现 三 种 方法 得 到 的 解 的 结构 大 致 相似 ， 与 参考 文献 中 的 结果 一 致 。 图 2 
给 出 了 三 种 方法 计算 结果 中 ， 在 y = 0 处 密度 和 总 磁场 强度 随 x 的 变化 。 从 中 可 
以 看 出 ， 在 激流 处 TVD-GLM-TD 格式 的 耗 散 最 小 ， 这 种 方法 有 效 的 降低 了 算法 的 
数值 粘性 。 
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图 1 二 维 Rotor 问题 的 密度 p 和 磁场 B? 在 t = 0.15 时 刻 的 等 值 线 .(aj 和 (b) 采 用 的 TVD-LF 方法 ， 


(cj) 和 (d) 采 用 TVD-GLM 修正 ，(e) 和 (f) 采 用 TVD-GLM-TD 修正 。 


Fig. 1 Level curves of density p and magnetic field magnitude B? of 2D Rotor 


question at t — 0.15. 
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图 2. —1 Rotor 问题 : t = 0.15 时 刻 x 方 向 ， 使 用 不 同方 法 计算 的 密度 p 和 磁场 |B| 对 比 图 。 


0.2 0.4 


Fig. 2 Comparison diagram of density p and magnetic field magnitude |B| by using different 


methods in x direction at t = 0.15. 


四 、 磁 云 边界 问题 数值 试验 
1、 计 算 初 始 条 件 

磁 云 是 行星 际 空间 中 磁场 和 等 离子 体 的 一 种 大 尺度 扰动 现象 。 其 特征 是 较 低 
的 等 离子 体 密度 和 温度 、 较 低 的 等 离子 体 B 和 磁场 方向 的 旋转 。 对 行星 际 磁 云 的 
研究 有 助 于 了 解 CME 爆发 时 物质 质量 、 磁 场 和 能 量 的 相互 关系 。Fan[13] 通 过 数 
值 模拟 磁 云 和 电流 片 相互 作用 过 程 ， 验 证 了 Wei[11] 对 磁 云 边界 层 给 出 的 定义 。 
为 了 用 空间 物理 研究 中 的 实际 问题 验证 改进 算法 的 计算 效果 ， 本 文 对 这 一 问题 ， 
用 [13] 中 的 TVD-LF 原始 方法 和 新 改进 的 GLM-TVD-LF 和 GLM-TVD-TD 两 种 方法 进 
行 了 数值 模拟 。 模 拟 采 用 和 矩形 网 格 ， 网 格 间距 Ax = 2 x Ay = 0.2， 计 算 区 域 


Lx X Ly =20x20。 令 x = —10.0 +i x Ax,y = —10.0 +j x Ay， 初 始 的 背景 磁场 为 


Harris 型 电流 片 : 


By = B,tanh (* Ax) i (17) 


磁 云 初始 的 中 心 位 置 (xu yo) = (0,3.5), RÁTRVEAR ER, = 5000， 磁 云 的 位 形 采 
用 Lundquist[20] 给 出 的 无 力 场 模型 : 


当 7 S Te 时 ， B, = 0, Bo EE CoBoJ1 (ar), B; = CoBoJo (ar); 


rp—r r2 rp 


sS o c-r 时， B, - B, (1— 


) cosg, Bo = —Bo (1 十 $) sing, 


7 一 re r? 7 一 rc 


B, = 0。 其 中 ，r 为 计算 点 距离 磁 云 中 心 的 距离 ， 磁 云 半径 xr, = 3.0，7 = 20r- 


Jo fJ 2&8 1 X NÆ Besse) KA a NENIAE Ha = 一 2.4/r, Cy = 5.7850. 


2、 数 值 模拟 结果 


首先 通过 数值 试验 确定 GLM-TVD-TD 修正 中 可 调 参数 = 的 最 佳 取 值 。 图 3 为 


计算 进行 到 t = 2 时 刻 , 不 同 耗 散 修 正 系数 e 模 拟 结果 的 磁场 强度 |B| 和 等 离子 体温 
度 T 的 模拟 结果 在 x = 0 处 随 y 变化 的 曲线 。 计算 结果 表明 当 e < 0.5 时 , 如 e = 0.3 
时 , 得 到 的 曲线 在 间断 面 和 极 值 附近 有 细微 的 振荡 出 现 , 而 e = 0.5 对 应 的 结果 过 


EKRE, 在 计算 稳定 的 同时 具有 较 低 的 数值 粘性 .因此 在 此 问题 的 数值 模拟 中 ， 


选取 e = 0.5. 


IB| 
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图 3 !=2 时 刻 ， 不 同 耗 散 修 ] 
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FE 系 数 e 对 应 磁场 |B| 曲 线 和 等 离子 体温 度 7 曲线 。 


Fig. 3 Curves of magnetic field magnitude |B| and plasma temperature T by using 


different diffusion reduction coefficie 


nt g at t= 0.15. 


图 4 给 出 了 三 种 方法 的 模拟 结果 在 上 = 2 时 刻 ， 沿 着 x = 0 方向 跨越 磁 云 边界 
层 时 磁场 强度 |B|、 磁 场 x 轴 夹 角 中 p、 等 离子 温度 T、 等 离子 体 密度 D 以 及 等 离子 
热 压 和 磁 压 之 比 的 对 数 logB 随 y 的 变化 曲线 。 对 照发 现 ， 三 种 方法 得 到 的 结果 图 
形 大 致 相同 。 但 是 相 比 于 [13] 中 TVD-LF 方法 得 到 的 结果 ，GLM-TVD-LF 和 


GLM-TVD-TD 方法 的 结果 中 y = 
TVD-LF 方法 求解 8 波形 式 的 M 


一 2.5 附 近 的 重 联 点 的 位 置 有 所 变化 。 我 们 在 应 用 
HD 方程 ， 再 次 对 本 问题 进行 计算 后 ， 发 现 8 波 法 


结果 中 重 联 点 的 位 置 与 GLM-TVD-LF 及 GLM-TVD-TD 的 结果 更 接近 。 文 献 [13] 使 用 
柏 松 校正 法 控制 磁场 散 度 ， 并 用 一 种 伪 时 间 步 方法 求解 柏 松 方程 。 我 们 的 数值 试 
验 表 明 这 种 求解 柏 松 方程 的 方法 并 不 能 完全 收敛 。 这 可 能 是 造成 重 联 点 位 置 变化 


的 原因 。 同 时 ， 在 磁 云 中 心 ，GLM-TVD-LF 和 GLM-TVD-TD 格式 的 磁场 方向 的 变 
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化 区 更 罕 ， 说 明 这 两 种 方法 的 数值 耗 散 较 低 。 
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图 4 T = 2 时 刻 ， 沿 着 x = 0 使 用 不 同方 法 模拟 磁场 强度 |B|、 磁 场 和 x 轴 夹 角 中 5 、 等 离子 温 
度 T、 等 离子 体 密度 D 以 及 等 离子 热 压 和 磁 压 之 比 的 对 数 logB 随 y 的 变化 曲线 。 图 中 实 线 使 用 
的 是 求解 GLM-MHD 的 TVD-TD 方法 ， 点 划 线 使 用 的 是 求解 GLM-MHD 的 TVD 方法 ， 虚 线 使 
用 的 是 [13] 中 的 原始 数值 方法 。 


Fig. 4 Variety curves of magnetic field magnitude |B|, the angle between the filed vector B and 
the x-axis Pg, plasma temperature T, plasma density D and the ratio B of the magnetic 


pressure and plasma pressure. 


为 了 定量 观察 磁场 散 度 V . 互 的 大 小 和 变化 情况 , 本 文 计 算 了 平均 磁场 散 度 误 
差 [21]: 


Xia), | 
M 


error — (18) 


AP M 为 总 的 网 格 数 。 


图 5 给 出 用 TVD-LF. TVD-GLM 和 TVD-GLM-TD 方法 对 磁 云 激 波 相 互 作 用 问题 
进行 模拟 时 ,计算 结果 的 平均 磁场 散 度 误差 随时 间 变 化 曲线 。 从 图 中 可 知 ,， 三 种 
方法 中 , TVD-LF 方法 使 用 泊 松 校正 法 ,对 磁场 散 度 误差 的 抑制 最 为 有 效 ， 其 次 是 
TVD-GLM 方法 。 随 着 计算 的 进行 ， 三 种 方法 得 到 的 平均 磁场 散 度 误差 均 逐 渐 减 
小 。 应 用 GLM 方法 控制 磁场 散 度 误差 时 ， 虽 然 磁场 散 度 误差 绝对 值 比 泊 松 校正 
法 要 大 , 但 前 文 展示 的 计算 结果 表明 磁场 散 度 误 差 对 计算 结果 无 明显 影响 。 这 说 
明 在 本 算 例 中 ，GLM 方法 能 够 有 效 的 控制 磁场 散 度 误差 。 


log(error) 


图 5 error 随时 间 变 化 曲线 


Fig.5 The variety curves of error over time 


本 算 例 的 数值 程序 使 用 Intel Fortran 编译 器 ， 以 -02 优化 选项 编译 ， 在 配备 
Intel Xeon E7450 CPU(2.4GHz 主 频 ) 的 计算 机 上 运行 。 运 行 到 物理 时 间 t=2 时 刻 ， 
文献 [13] 的 TVD-LF 格式 耗费 的 计算 时 间 为 307 秒 ， 而 本 文 的 GLM-TVD-LF 和 
GLM-TVD-TD 方法 中 耗费 的 计算 时 间 为 123 秒 。 文 献 [13] 的 TVD-LF 在 用 柏 松 校正 
法 控制 磁场 散 度 时 ， 需 要 用 迭代 的 方法 求解 柏 松 方程 ， 需 要 额外 耗费 较 多 时 间 。 
这 可 能 是 其 计算 耗 时 较 长 的 原因 。 

五 、 结 论 

本 文 对 2.5 维 笛 卡 尔 坐标 系 下 的 MHD 方程 组 采用 TVD 格式 差分 ， 分 别 对 耗 
艇 项 和 磁场 散 度 约束 条 件 进行 修正 ， 发 展 出 两 种 快捷 上 共有 TVD 特性 的 组 合 数值 
模拟 新 方法 。 有 别 于 其 它 的 TVD 格式 ，TD 修正 稳定 有 效 ， 可 对 耗 散 项 进行 人 为 
优化 ，TVD-GLM 方法 则 简单 易 行 ， 相 同 条 件 下 模拟 计算 耗 时 减少 一 半 。 在 二 维 


Rotor 问题 中 所 得 结果 与 相关 文献 吻合 ， 能 够 准确 模拟 复杂 螺旋 阿尔 芬 波 结构 。 
在 磁 云 和 电流 片 相 互 作 用 过 程 的 数值 模拟 试验 中 , 两 种 修正 方法 得 到 的 结果 具有 
更 高 的 精确 度 和 分 辨 率 ， 与 观测 数据 相符 。 为 今后 将 其 推广 到 更 为 复杂 的 高 维 数 
值 模拟 计算 莫 定 了 坚实 的 基础 。 
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